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Abstract: We construct a new gravitational action which includes cubic curvature in- 

m ■ 

\f} • teractions and which provides a useful toy model for the holographic study of a three 

£T) . parameter family of four- and higher-dimensional CFT's. We also investigate the black 

hole solutions of this new gravity theory. Further we examine the equations of motion of 
quasi-topological gravity. While the full equations in a general background are fourth-order 
in derivatives, we show that the linearized equations describing gravitons propagating in 
^ ; the AdS vacua match precisely the second-order equations of Einstein gravity. 
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1. Introduction 

Recently, there has been some interest in gravitational actions with higher curvature actions 
in the context of the AdS/CFT correspondence. For example, Einstein gravity in the AdS 
bulk defines a universality class of CFT's in which the ratio of the shear viscosity to entropy 
density is given by precisely rj/s = l/(47r) [1, 2]. However, it is understood that adding 
higher curvature interactions to the bulk gravity action leads to a broader class of CFT's 
in which this ratio generally depends on the value of the additional gravitational couplings 
[3]. Further it is possible with some holographic constructions to violate the famous bound 
conjectured by Kovtun, Son and Starinets (KSS) [1] producing theories with -q/s < l/(47r). 
In certain string theory constructions, the appearance of curvature-squared interactions 
produces violations of the KSS bound [4] but these models only produce reliable results in 
a regime where the corresponding gravitational coupling is parametrically small. Hence at 
present, one can only deviate perturbatively away from the universality class defined by 
Einstein gravity in these string theory models. 
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However, it is also of interest to explore situations where the gravitational couplings are 
finite. For example, holography can yield new consistency conditions for the gravitational 
theories and their dual CFT's. One theory which provides a useful toy model in this 
regard is Gauss-Bonnet (GB) gravity. Even with a finite coupling for the curvature-squared 
interaction, this theory still provides some calculation control, which has been exploited in 
several recent holographic studies [5, 6, 7, 8, 9]. However, GB gravity only introduces a 
single new coupling which limits the range of dual CFT's which can be studied. A natural 
generalization would be the further addition of interactions cubic in the curvature, as this 
allows the investigation of the the full range of parameters in the three-point function of 
the stress tensor [10]. A straightforward extension of GB gravity would be to include the 
cubic interaction of Lovelock gravity [14]. However, because of the topological origin of the 
Lovelock terms the cubic interaction only contributes to the equations of motion when the 
bulk dimension is seven or greater. In the context of the AdS/CFT correspondence, this 
means that such a term will be effective in expanding the class of dual CFT's in six or more 
dimensions [8, 9]. Our key result in this paper is to construct a new gravitational action 
with cubic curvature interactions which provides a useful toy model to study a broader 
class of four (and higher) dimensional CFT's, involving three independent parameters. In 
the following, we describe the construction for the new gravitational action and investigate 
black hole solutions in this theory. We leave the detailed study of the properties of the 
dual class of CFT's to a companion paper [10]. 

An outline of the rest paper is as follows: We begin with a review of black hole solutions 
and various aspects of these solutions in Gauss-Bonnet (GB) gravity coupled to a negative 
cosmological constant in section 2. Inspired by the GB equations of motion determining 
black hole solutions, we construct a new interaction that is cubic in curvatures and yields 
similar simple solutions in section 3. Again we wish to emphasize that this interaction is 
not the six-dimensional Euler density as appears in third-order Lovelock gravity. Further, 
we show in appendix A that the new interaction does not have a topological origin and 
hence we call the new theory: 'quasi-topological gravity' We turn to a discussion of the 
asymptotically AdS black hole solutions in section 4. While the focus of this discussion is 
planar black holes in five dimensions, we generalize the results to curved horizons and higher 
dimensions in sections 4.2 and 4.3. In section 5, we examine black hole thermodynamics 
in the new theory, deriving some of the basic thermal properties of the black holes and 
the corresponding plasmas in the dual CFT. We examine the equations of motion of quasi- 
topological gravity in section 6. While the full equations in a general background are fourth- 
order in derivatives, we show that the linearized equations describing gravitons propagating 
in the AdS vacuum solutions are precisely the second-order equations of Einstein gravity. 
We conclude with a brief discussion of our results and future directions in section 7. 

While we were in the final stages of preparing this paper, two related preprints appeared 
in which exceptional new theories of curvature-cubed gravity were constructed. Ref. [11] 
constructs an interesting curvature-cubed theory in three dimensions. Up to a contri- 
bution proportional to the six-dimensional Euler density, the curvature-cubed interaction 
constructed in five dimensions by [12] is identical to that studied here. Refs. [12, 13] are 
also able to relate our interactions in D > 7 to Weyl-invariant combinations of curvatures 
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combined with the six-dimensional Euler density. 



2. Black Holes in Gauss-Bonnet gravity 

We begin here with a brief review of black holes in Gauss-Bonnet (GB) gravity. The latter 
corresponds to a theory of gravity in which a curvature-squared interaction is added with 
the form of the density for the Euler characteristic of four-dimensional manifolds, 

<^4 = RabcdR ahCd — 4i? a fe-R a6 + B? ■ (2.1) 

Of course, this term will not affect the gravitational equations of motion if the dimension 
of the spacetime is four (or lower), however, it makes interesting contributions for D > 5. 
GB gravity can be seen as the simplest example of the Lovelock theories [14] discussed 
above. As explained, despite having a higher curvature action, the resulting equations 
of motion are still second-order in (time) derivatives and this produced some interest in 
early discussions of higher curvature corrections to string theory [15, 16]. These discussions 
also lead to an extensive study of black hole solutions in this theory [17]. More recently, 
there has been some renewed interest in asymptotically AdS black hole solutions in GB 
gravity [18, 19], especially in the context of the AdS/CFT correspondence [5, 6, 7]. In our 
following, we discuss the black hole solutions focussing on GB gravity with D = 5 and with 
a negative cosmological constant: 



I = — — — / d 5 x J^g 
16vrG 5 y V 



2 
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(2.2) 



We add some comments about higher dimensions at the end of this section. 

Let us present the ansatz for the metric of five-dimensional planar AdS black holes, 
which we will be using throughout the paper: 

ds 2 = -= (-N(r) 2 f(r) dt 2 + dx 2 + dy 2 + dz 2 ) + -_- dr 2 . (2.3) 
L A r j(r) 

Inserting this metric ansatz into the action (2.2) (and integrating by parts a number of 
times) yields 

where the 'prime' indicates differentiation with respect to r. Schematically, the equation 
of motion coming from varying the lapse N takes the simple form [r 4 (- ••)]' = and so the 
metric function / is given by solving for the roots of a quadratic polynomial [18, 19]: 

Xf{rf - f(r) + 1-^ = 0. (2.5) 

The latter yields two solutions 

f±(r) = ^ L± 1-A( I -~ T ! (2.0) 
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Now varying 5f yields a constraint which requires that N = constant, which we leave 
unspecified for the moment. In the following, we will consider only the solutions with /_, 
since the other branch with /+ contains ghosts and is unstable [16] — as we will see in 
later sections. With the choice / = /_, it is easy to verify that the horizon appears at 
r = r h =uj. 

In fixing the value of N, it is convenient to consider the solution with u = 0, 

r 2 I? 
ds 2 = -r (-iV(r) 2 /oo dt 2 + dx 2 + dy 2 + dz 2 ) + dr 2 . (2.7) 



where we have adopted the notation 



1 



/oo = lim f(r) = — 1 - VT^iX • (2.8) 

r— >oo Za L 

We recognize eq. (2.7) as anti-de Sitter (AdS) space, presented in the Poincare coordinates. 
From g rr above, we also see that the AdS curvature scale is given by L = L/y/J^. This 
metric also makes apparent a convenient choice for the lapse, namely A 2 = l//<x>, which we 
adopt in the following. In the AdS vacuum (2.7), this ensures that any motions in the brane 
directions are limited to lie within the standard light cone, i.e., = — dt 2 + dx 2 + dy 2 + dz 2 . 
In the black brane solution (2.3), we still have lim r ^oo A 2 f(r) = 1 and so this comment 
applies in the asymptotic region. In the context of the AdS/CFT correspondence, the 
latter means that the speed of light in the boundary CFT is simply c = 1 . 

Examining the solutions in eq. (2.6) or in eq. (2.8), we see that there is an upper 
bound at A = 1/4. For larger values of A, the gravitational theory does not have an anti-de 
Sitter vacuum and the interpretation of the solutions (2.6) becomes problematic. In fact, 
using the AdS/CFT correspondence to demand consistency of the dual CFT, e.g., requiring 
that the boundary theory is causal, imposes much more stringent constraints on the GB 
coupling [5, 6] 

7 9 

< A < . (2.9) 

36 - - 100 K ' 

We now turn to the thermodynamic properties of these GB black holes. The temper- 
ature of the black brane solutions is then given by the simple expression: 



1 rlf 



4vr L 2 ttL 2 
which when evaluated with N = l/^/f^ becomes: 



rh N = ^N, (2.10) 



7T 



L 2 



1 ( 1 + VT^4A) 



1/2 

(2.11) 



The latter can be calculated by the standard technique of by analytically continuing the 
metric (2.3) to Euclidean time, r = —it, and choosing the periodicity of r to ensure the 
geometry is smooth at = uj. Next we evaluate the Euclidean action: 
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where V3 is the regulator volume obtained by integrating the (x, y, z) directions. Further 
we have limited the radial integration from r = uj to r + to regulate the asymptotic or UV 
divergence in Ie- The divergent r\ contribution is removed with background subtraction 
using the AdS vacuum (2.7), i.e., that is Ie — 1% remains finite in the limit r + — > 00. Note 
that in general, such a calculation would include a generalized Gibbons-Hawking term [20], 
as well as other boundary terms to regulate the divergences in the Euclidean action [21]. 
However, these surface terms do not contribute to the final result for planar AdS black 
holes when we use the background subtraction approach. Therefore we identify the free 
energy density as 

Then we may identify the energy and entropy densities as 



P--T^/T)= - — / , (2.14) 



UJ 3 



dT 4G 5 4G 5 L 3 ' 

One can confirm that the last result matches the entropy calculated using Wald's techniques 
[22]. 

These solutions are easily generalized from five to an arbitrary spacetime dimension, 
D. In this case, the action is conveniently parameterized as 



(D- 1)(D- 2) , XL 2 

T 2 +R+ (D-3)(D-4f\ 



(2.15) 



We also generalize the metric ansatz to include spherical and hyperbolic, as well as planar, 
horizons: 

v 2 civ 2 
ds 2 = -(k+ f( r ))N(r) 2 dt 2 + + r 2 dl\ (2.16) 

L k + j^fir) 

where dl\ is given by 

k = +1 : dn 2 D _ 2 (metric on S D ~ 2 ) , 

k = 0: tVE(^) 2 > (2-17) 

i=i 

k = —1 : dY? D _ 2 (metric on H D ~ 2 ) . 



Note that for k = ±1, the above line element has unit curvature. With this general ansatz 
incorporating both curved horizons and D spacetime dimensions, / is determined by simply 
solving for the roots of 



D-l 



A/(r) 2_ /(r) + 1 _^__ = 0) (2 . 18) 

and the solutions take the form given in eq. (2.6) with the replacement w 4 /r 4 — > lj®^ 1 /r ^ 1 . 
The lapse is again a constant and as above we choose N = l/V/oo- in general, the hori- 
zon is determined by fWh) = —k^ and so we only have = uj for the planar horizons, 
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i.e., k = 0. In the case of the curved horizons, explicitly evaluating requires solving a 
(D — l)-order polynomial in r^. Hence we only have a relatively simple solution in five 
dimensions where: 




(2.19) 



3. New Curvature-Cubed Interaction 

As discussed in the introduction, we are motivated by considerations of the AdS/CFT 
correspondence to consider a curvature-cubed theory of gravity in five dimensions. GB 
gravity has a number of features which one might want to reproduce, such as providing 
second-order equations of motion and a family of exact black hole solutions. A natural 
candidate to extend these properties to a curvature-cubed theory would be the Lovelock 
theory where the six-dimensional Euler density is added as a new gravitational interac- 
tion. However, this curvature-cubed interaction would only contribute to the equations of 
motion in seven and higher dimensions and hence will not contribute in the desired five 
dimensions. While Lovelock's work then indicates that it should not be possible to find 
an alternate action which yields second-order equations of motion, we begin by writing 
the most general interaction including all possible curvature-cubed (or more precisely, six- 
derivative) interactions in five dimensions and attempt to tune the coefficients to produce 
a simple equation for the black hole solutions, as discussed for GB gravity in the previous 
section. We return to the equations of motion in section 6 and we demonstrate that the 
linearized equations of motion in the AdS vacuum are indeed second-order. 
Let us begin by listing a basis of the possible six-derivative interactions: 



In assembling this list, we have discarded any total derivatives, e.g., X7 a X7 a X7 b X7 c R bc and we 
have simplified various expressions using the index symmetries of the Ricci and Riemann 
tensors. In particular, these symmetries allow us to reduce any other index contraction of 
three Riemann tensors to some combination of terms 1 and 2. Further, term 12 can be 
reduced to term 13 using V a R ab = ^V b R. Similarly, using the Bianchi identities, terms 9 
and 10 can be shown to be reducible to other terms and total derivatives as well. Hence, we 
are left with a list of 10 independent interactions which are cubic in curvatures. Combining 
all of these interactions together in a single expression gives: 



q p r>abc pi 
A P nabcd p 

5. R abcd R ac R M 




6. R a b R b c R c a 

7. R a b R b a R 

8. R 3 

9. V a R hcde V a R bcde 

10. V a V c R abcd R bd 



11. V a R bc V a R> 

12. V a R ab V b R 

13. V a RV a R . 



- V~9 (ci Rab dR cd R ef b + c 2 Rab^Rcd^ R ef ab + c 3 RabcdR^eR^ 

+ c 4 R a bcdR abcd R + c 5 R abc dR ac R bd + c 6 R b R b R^ + c 7 R b R b R 




(3.1) 
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At this point, we substitute the black brane metric ansatz (2.3) and evaluate eq. (3.1). 
The next step will be to see if the coefficients c { can be tuned to produce a result with the 
same form as in eq. (2.6). In order to accomplish this task, we integrate by parts repeatedly 
to put as many terms as possible in the form N(r) x (• • • ) where the factor in brackets is 
independent of the lapse function. The resulting expression then becomes 



—7 N 

-gZ -- 



(8r 9 (c„ + 2c 13 ) / 2 / {6) + (r 8 (236c n + 496c 13 ) f 2 + r 9 (56c 13 + 28c 13 ) 



4L 9 

x //')) / (5) + {r 9 (24c 3 + 48c 2 + 20c n + 48c 8 + 12c 6 + 24c 7 + 12c 5 + 48c 4 
+c 13 ) //" + 8r 9 (c n + 2c 13 ) ( f) 2 + r 7 (5056c 13 + 384c 4 + 96c 6 + 2264c n 
+96c 2 + 96c 3 + 288c 7 + 960c 8 + 72c 5 ) f 2 + r 8 (72c 5 + 288c 4 + 120c 3 
+1312c 13 + 192c 2 + 480c 8 + 192c 7 + 608c„ + 84c 6 ) //') / (4) + (r 8 (48c 4 
+12c 5 + 48c 8 + 24c 7 + 16c 13 + 12c 6 + 8c n + 24c 3 + 48c 2 ) / {f"f 
+ (V 8 (240c 8 + 96c 7 + 60c 3 + 100c n + 96c 2 + 224c 13 + 36c 5 + 144c 4 + 42c 6 ) (f) 2 

+r 7 (1494c 5 + 6069c 4 + 3054c 2 + 4380c 7 + 4012c n + 9776c 13 + 12000c 8 
+1794c 6 + 2340c 3 + 18d) //' + (r 9 (6c 5 + 24c 2 + 24c 4 + 12c 7 + 24c 8 + 12c 3 
+6c 6 ) /' + r 8 (1296c 2 + 920c 13 + 384c B + 708c 3 + 388c„ + 2016c 8 + 1536c 4 
+414c 6 + 888c 7 ) /) /" + r 6 (1584c 6 + 8140c n + 36c x + 1440c 2 + 1212c 5 
+4944c 7 + 1488c 3 + 5952c 4 + 17280c 8 + 19712c 13 ) f 2 ) f" - 2r 9 (c B + 4c 8 
+c 6 + 2c 7 + 4c 4 + 2c 3 + 4c 2 ) {f'f + (r 8 (42c 3 + 84c 8 + 21c 6 + 84c 2 + 42c 7 
+21c 5 + 84c 4 ) /' + r 7 (10896c 8 + 4260c 7 + 18d + 6384c 4 + 1308c n + 3696c 13 
+1608c 5 + 2700c 3 + 4488c 2 + 1842c 6 ) /) (f"f + (r 6 (26544c 13 + 306d 
+28092c 7 + 8790c 5 + 84000c 8 + 10740c 6 + 8916c n + 17472c 2 + 12600c 3 
+35088c 4 ) //' + r 5 (324 Cl + 5436c 5 + 564c 2 + 6960c 6 + 82560c 8 + 8676c n 
+22608c 7 + 6069c 3 + 24384c 4 + 24192c 13 ) f 2 + r 7 (264c n + 252c B + 1056c 4 
+720c 7 + 300c 6 + 624c 2 + 1920c 8 + 672c 13 + 408c 3 ) (/') 2 ) /" - 16r 3 (8c 3 

+40c 4 + 400c 8 + 80c 7 + 16c 6 + 3d + 16c 5 + 4c 2 ) f + r 6 (434c 3 

+1240c 4 + 301c 5 + 950d + 592c 2 + 2800c 8 + 9d + 361c 6 ) (/') 3 + r 5 (10122c 6 

+10236c 3 + 558d + 10080c 13 + 2160di + 29040c 7 + 12480c 2 + 8202c 5 

+31296c 4 + 96000) (ff + r 4 (-900di + 5904c 3 + 24240c 7 + 504 Cl 

+23520c 4 + 7248c 6 + 5748c 5 + 91200c 8 + 5232c 2 )) f 2 f] +■■■. (3.2) 

Note that not all terms can be put in the desired form with further integration by parts 
and so the '• • • ' indicates the presence of spurious terms containing factors like (N") 2 /N, 
for example. Focusing on the terms appearing explicitly in eq. (3.2), we find that choosing 
the values of the q's as 

9 60 33 54 

i- c 3 — —j c i y ° 2 °7 — — Ci y ° 2 
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3 3 15 11 .„„, 

2. c 4 = -d + -c 2 6. c 8 = — d + — c 2 (3.3) 

15 72 

3. C 5 =yC 1 + yC 2 7. dl = 

18 64 

4- C 6 =yC 1 + yC 2 8. d 3 = 

reduces this expression to the following simple form 

v _ z= l|^) (ci + 2c2)(r4/3)/ _ (34) 

At the same time, the spurious terms denoted by '• • • ' in eq. (3.2) also vanish with this 
choice of coefficients. It is quite remarkable that there is enough freedom in the general 
action (3.1) to produce this simple result. In fact, we are still free to choose the (relative) 
values of d and c 2 in constructing this interaction. Explicitly then, if we choose c- L = 
1, c 2 = 0, the new curvature-cubed interaction takes the form 



pcdpe/pab, ^ /ni n r>abcd r> 70 E> r>abc ride 
K ab K cd K ef +-^[ ZLK abcdtt K ~ <2 H abcd K e K 



+120 R abcd R ac R bd + 144 R a b R h c R c a - 132 R b R b a R + 15 i? 3 ) (3.5) 
or with d = 0, c 2 = 1, 

Z' 5 = R ab cd R c /fR ef ab + ^ (21 R abcd R abcd R- 120 R abcd R abc e R de 

+144 R abcd R ac R bd + 128 R a b R b c R c a - 108 R b R b a R + 11 i? 3 ) . (3.6) 

The fact that we do not produce a unique interaction should not be surprising. Any 
curvature-cubed interaction can be modified by the addition of the six-dimensional Euler 
density X§ without affecting the equations of motion. In fact, we can infer the form of X% 
by setting c x = — 2c 2 , in which case eq. (3.4) vanishes, as it must if evaluated for the six- 
dimensional Euler density. A standard normalization for the six-dimensional Euler density 
is: 

<*6 = g Zabcdef E 9hljkl Rab^ Rcd %j Ref H 

a Tf cdjj e f tj ab QTjcdr/e.fjyab n/lf? r>abc jjde , o r> r>abcd p 

— 4/t ab K cd K e f —8K ab K cd K e f — l^ti abcd ti e K + 6 K abcd K ti 

+24 R abcd R ac R bd + 16 R b R b R c a - 12 R b R b a R + R 3 , (3.7) 

where in the first line, e abcde f is the completely antisymmetric tensor in six dimensions and 
hence the corresponding expression only applies for D = 6. However, the first line also 
makes clear that this expression should vanish when evaluated in five (or lower) dimensions. 
This normalization corresponds to the choice c 2 = 4 and d = —8. That is, X% = 4Z' 5 — 8Z^,. 

3.1 Generalizing to D > 5 

At this point, we turn to generalizing this construction to higher dimensions. Given the 
freedom discussed above, we begin by setting d = 1 and c 2 = 0. Comparing to eq. (3.7) 
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shows we are guaranteed that, if a nontrivial interaction exists, the result will be distinct 
from the six-dimensional Euler density. With this choice, we substitute the L>-dimensional 
extension of eq. (2.3) into the action (3.1). One then finds with a judicious choice of the 
remaining q's, the result can be reduced to 

V^Z D ~ ^ (r D -'f{rf)' . (3.8) 

The required choice of coefficients (with d = 1 and c 2 = 0) is: 

1. c,(D) = - 3(^-2) ,,*<*>--; 3(3I> - 4) 



(2D-3)(D-4) ' v 7 2(2L> - 3)(D - 4) 

2 C * (D) = 8(20-V-4) C «<°> " 8(2D-fl(fl-4) < 3 ' 9 > 

3D 

3 C - (P)= (2 D - 3 )( D -4) 7. c„( D ) = 

4 *< P >= (2P-V-4) *<*(«) = <>. 

In practice, we determined the coefficients separately for D = 5 ... 10 and found the general 
expressions above to fit the results in all of these cases. Given these expressions, the general 
form of Zr> becomes 

-7 —Rcdpe/pabj 1 ( 3(3ZJ - 8) bcd 

Zu-R ab R cd R ef + {2D _ 3){D _ 4) [ ~ 8 R abcdR R 

- 3(D - 2)R abcd R abc e R de + 3DR abcd R ac R bd (3.10) 
+ 6(D - 2)R a b R b *R c ° - 3(3D ~ 4) R a b R b a R + ^R 3 

One easily verifies that this result reduces to eq. (3.5) for D = 5. 

An important note is that, with the coefficients prescribed above for D = 6, the 
resulting action (3.8) is trivial, i.e., there is an overall factor of zero in this result. Hence 
Zq does not actually produce a nontrivial interaction which is cubic in curvatures. One 
might be tempted to believe instead that Zq yields another topological invariant in six 
dimensions. Of course, there is no obvious known invariant to which Zq might correspond 
[23]. In appendix A, we demonstrate that f S'x^fgZQ is not a topological invariant by 
explicitly evaluating this expression for some nontrivial six-dimensional geometries. Hence, 
we refer to the theory of gravity extended with our new curvature-cubed interaction as 
'quasi-topological gravity' At this point, we also note that our construction does not yield 
a nontrivial curvature-cubed interaction for D < 4. 

Of course, one can generalize the interaction (3.10) for D > 6 by adding another com- 
ponent proportional to the six-dimensional Euler character (3.7). This would be equivalent 
to leaving c 2 arbitrary in our analysis above. Hence, we complete the discussion by gen- 
eralizing eq. (3.6) to higher dimensions using the formula: Z^ = 2Z D + \Xq. The final 
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expression can be written as 



Z' D = R ab cd R cd e fR ef ab + {2D _z ){D _ A) (" 12 (° 2 -5^ + 5) R abcd R abc e R de 

+ l (£>2 _ 4jD + 2 ) RR abcd R abcd + 12 (D - 2) (D - 3) R abcd R ac R bd (3.11) 
+8 (.D - 1) (D - 3) R a b R b c Rc a - 6 (Z> - 2) 2 i?i? a % a + \ (D 2 - AD + 6) R 3 ) . 



4. Black Hole Solutions 

We have thus far constructed a new gravitational action that includes interactions up to 
cubic order in the curvature and which still yields particularly simple equations to find 
black hole solutions. In this section, we complete the study of the black holes in this new 
theory. While we have written an action for the theory in arbitrary number of spacetime 
dimensions, we will focus on the case D = 5 here. Further we begin by examining solutions 
of the form given in eq. (2.3) and hence construct black holes with planar horizons [24, 18]. 
The extension of this analysis to larger D and curved horizons (with spherical or hyperbolic 
geometries) is straightforward and will be discussed briefly at the end of this section. 

4.1 Planar Black Holes 

We begin with the five-dimensional action: 



16vrG 5 J 



12 XL 2 7uL A 



— + R+—X 4 + ^ r Z 5 



(4.1) 



which extends the GB action with the addition of the curvature-cubed interaction Z§. 
Next, as in section 2, we consider the following metric ansatz: 

r 2 1? 
ds 2 = (-N(r) 2 f(r) dt 2 + dx 2 + dy 2 + dz 2 ) + -^-^y dr 2 . (4.2) 

Evaluating the action (4.1) with this metric then yields 

where the prime again denotes a derivative with respect to r. The variation 5N now yields 

[r 4 (l-/ + A/ 2 + ^ 3 )]' = 



l-/ + A/ 2 + / u/ 3 = ^. (4.4) 



UJ 



4 



r 



Similarly, satisfying the equation produced by taking the variation of / requires that either 
N' = or -1 + 2A/ + 3/// 2 = 0. Since the latter is generally inconsistent with eq. (4.4), we 
arrive at N = constant. As in section 2, we choose N 2 = 1/foo where foo = lim,.-^ f(r). 
This choice ensures that the speed of light in the boundary metric is just one. 
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We are now left with a cubic equation (4.4) to solve for f(r). To do so, we first make 
the substitution / — x , with which eq. (4.4) becomes: 

^-»(^). + »( ^±g^ )-o. («, 

This expression is further simplified by defining 



3/x + A 2 2A 3 + 9^ + 27^(1-^) 
^=^2- Q = 54^ • (46) 



We then arrive at the depressed form of the equation: 

x 3 - 3px - 2q = 0. (4.7) 

In the following discussion, note that the r-dependence is entirely contained in the coeffi- 
cient q. Before proceeding, we observe that there are three distinct cases depending on the 
sign of the discriminant, V = q 2 — p 3 , of eq. (4.7) with the following results: 

1. q 2 — p 3 > =>• 1 real root and 2 complex roots conjugate to one another 

2. q 2 — p 3 < 3 unequal real roots 

3. q 2 — p 3 = =>■ 3 real roots, at least 2 of which must be equal 
Assuming that p ^ 0, we define 

a=(g+ V ^^) 1 , (4.8) 

i 

P = (q- VQ 2 ~ P 3 ) 3 , 

which allows the roots of eq. (4.7) to be written in the simple form using Cardano's formula. 
Shifting these roots as above yields the following solutions: 

h = <* + 0--£-, 
6ijl 

h = -r 2 (a + P) + i%c,-P)-^, (4.9) 

A = ~(<» + fl-^<«-ffl-£. 

If V = q 2 - p 3 > 0, a and ft can be taken as real and x\ corresponds to the single real 
root. In this regime, the solution is then given by / = x\ — If V = q 2 — p 3 < 0, a and 
P are necessarily complex but implicitly eq. (4.9) still yields three unequal real roots. In 
this case, a and ft can be chosen to have conjugate phases, i.e., 

a = jp~ e i9/3 and p = y/p e' id/3 (4.10) 
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where cos# = q/p? and sin# = \/p 3 — q 2 /p2 . Here, we are using the fact that p is always 
positive in this domain. The solutions may then be cast in the explicitly real but implicit 
form: 

A = 2^008---, 



/ 2 = -v^(cos^ + V3sin0 - A (4.11) 
h = -y/p ^cos ^ - a/3 sin Q 



While the precise form of /(r) is determined by eqs. (4.9) and (4.11), these results offer 
little insight into the physical properties of the corresponding solutions, e.g., which solutions 
actually correspond to black holes. However, we will see below that much of the physics 
can be inferred directly from the cubic equation (4.4). 

At this point, it is convenient to consider the AdS vacuum solutions. As discussed in 
section 2, the latter can be found by setting f(r) to be constant (i.e., setting uj = 0) or 
alternatively taking the limit r — > oo. Setting to = and f(r) = foo in eq. (4.4) yields: 

Hfoo) = 1 - foo + Xfl + Hfl, = • (4.12) 
With the choice N 2 = l//oo, the five-dimensional metric (2.3) becomes 

ds 2 = (-dt 2 + dx 2 + dy 2 + dz 2 ) + , (4.13) 

which corresponds to the metric for AdS 5 in Poincare coordinates. Further from the g rr 
component, we see that the radius of curvature of the AdS 5 spacetime is 

L 2 = L 2 //oo • (4.14) 

Implicitly, we have assumed > 0, which is not always the case — see discussion below. 

Now let us examine these solutions as functions of the couplings, i.e., in the [i — A 
plane. In particular, if we insert the asymptotic limit r — > 00, the discriminant of eq. (4.7) 
is useful in determining the number of vacuum solutions at different points in the parameter 
space. As seen in eq. (4.6), p is unchanged in this limit since it is independent of r but q 
is slightly simplified: 

3fi + X 2 2A 3 + 9//A + 27^ 
P= ^~ q ^ = 54^ • (415) 

Given the previous discussion, we see that eq. (4.12) yields 3 real solutions for = 
q 2 ^ — p 3 < and 1 real solution for > 0. The vanishing of the discriminant, Poo = 0, 
reduces to a quadratic equation for fi with solutions 

" - h ~ I ± I " - 3A > 3/2 ■ < 416 > 

Eq. (4.16) generates the two (upper) curves in the //-A plane shown in figures 1 and 2. In 
the region bounded by these two curves, < and there are three vacuum solutions 
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Figure 1: The red and blue curves indicate the positive and negative branches of eq. (4.16), 
respectively, where T> — 0. The region bounded by these two curves is where T> < 0. The green 
curve indicates p = 0. The letter labels refer to the various regions described in table 1. The blue 
shaded region indicates those couplings for which there exist asymptotically AdS black holes. 

while outside of these two curves, > with one vacuum solution — as long as n ^ 0. Of 
course, /i = is a special axis in the parameter space corresponding to the Gauss-Bonnet 
theory discussed in section 2. Note that the positive branch of eq. (4.16) crosses the A-axis 
at A = 1/4, which is precisely the critical coupling in the GB gravity. Recall that in GB 
gravity for A < j, there are two vacua while no vacuum solutions exist for A > \ . 

Below we will find that another interesting boundary in the fi-X plane is p = 0, which 
corresponds to the lowermost (green) curve in figures 1 and 2 - note that p = is always 
below both branches of T> = 0, except at the points (A, fi) = (0,0) and (§,— From 
eq. (4.15), we see that p = simply corresponds to = — 4p. A distinguishing feature of 
p = is that eq. (4.5) becomes a perfect cubic equation. 

As a point of clarification, we should note that when eq. (4.12) yields real roots, the 
value of foo may be either positive (as assumed above) or negative. For example, it is easy 
to see that if fj, > one of the roots must be negative (since /i(/oo — > — oo) ~ fif^ < 
while /i(/oo = 0) = 1 > 0). Consistency demands A^ 2 > and so we can not use the same 
choice for the lapse when < 0. If instead we choose A^ 2 = l/|/oo|, the metric (2.3) 
becomes 

ds 2 = -2 (dt 2 + dx 2 + dy 2 + dz 2 ) (4.17) 

which corresponds to a particular set of coordinates on five-dimensional de Sitter space, 
where we observe that r plays the role of time. The radius of curvature of this dS 5 spacetime 
is L 2 = L 2 /!/^!. Even though with < 0, eqs. (4.9) or (4.11) still yield nontrivial 
solutions f(r), however, the corresponding metrics should be interpreted as (singular) 
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Figure 2: A closer examination of the lower right quadrant of figure 1, covered by the regions 
denoted (c) and (g) in table 1. All three curves intersect at (A,/i) = (§, — ^j), where the two 
branches of eq. (4.16) end. Region (g) is bounded by these two branches. We divide (c) into: (cq) 
where p < 0, (ci) and (C2) with p > and to the left and right of (g), respectively. 

cosmological solutions rather than black holes. Hence we will not consider solutions with 
/oo < further here. 

We summarize our results in table 1, which enumerates the various kinds of vacuum 
solutions in different regimes of the parameter space. Note that in this table, we have 
categorized the AdS vacua as either stable or 'ghosty,' i.e., whether or not the graviton is 
a ghost in a particular AdS vacuum. Recall that for GB gravity, we discarded one branch 
of the solutions in eq. (2.6) because the analysis of [16] showed that the graviton was a 
ghost in these backgrounds. A similar analysis applies for quasi-top ological gravity, as we 
will see in section 6. The key factor distinguishing this feature of the various vacua is the 
slope of the cubic equation determining : 

Woo) = -l + 2A/oo + 3/i/£. (4.18) 

In section 6, we show that this expression appears as a pre-factor in the kinetic term for 
gravitons propagating in a given AdS vacuum. The kinetic term has the usual sign when 
h'(foo) < and the wrong sign when /i'(/oo) > 0. Hence the stable or ghost-free AdS vacua 
in table 1 are distinguished by having h'(foo) > 0. Given that this factor is simply the 
slope of h(foa), it is easy to see that since /i(/oo = 0) = 1 then if there is one AdS vacuum 
(i.e., one root with > 0), it will be ghost-free. Similarly if there is more than one AdS 
vacuum, one of these will contain ghosts. 

To gain some further insight into the black hole solutions, we return to eq. (4.4) which 
we re-write as 

Hf) = (l " ^4 ) " / + A/ 2 + iif = • (4.19) 
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Table 1: Table of various vacua and black hole solutions. The column labeled 'BH solution' 
indicates which root (4.9) yields the nonsingular black hole solution. In case (c), the black hole 
solution is only realized in the regions denoted (c ) and (ci) in figure 2. 



Now in the asymptotic limit r — > oo, we simply recover eq. (4.12) and the roots match 
the vacuum solutions /oo- We can regard the effect of r decreasing through finite values 
as reducing the 'constant' term in the cubic polynomial of / and as a result, the roots of 
h(f) shift away from /oo- 

To illustrate various possibilities, figure 3 plots an example of h(f) in case (g) with 
three AdS vacua - see table 1. First we consider the smallest root ji, which asymptotically 
reaches the AdS vacuum solution with the smallest value of /oo- As shown when r decreases, 
this root decreases moving monotonically to the left until it reaches fi = at r = to. As 
r shrinks to even smaller values, / becomes negative and the solution becomes singular 
with /2 — > — oo as r — > 0. Of course, this behaviour is precisely that of a black hole with a 
horizon at r = u. 

Next consider the second root in figure 3. In this case, the slope h'(f) is positive 
and so the corresponding AdS vacuum contains ghosts. Note that as r decreases (or w 4 /r 4 
increases), the root now moves to the right, i.e., grows as we move to the interior of 
the solution. This behaviour is problematic as it corresponds to a negative mass solution 
and it seems to be connected to the ghost problems. As we discuss below, the solution 
reaches a naked singularity at r = ro (= 1.150a;, in this particular example) where the 
two roots, /i and fa, coalesce. One could overcome the problem with negative masses by 
simply choosing the integration constant oj 4 < 0. In this case, f% moves to the left with 
decreasing r but a naked singularity is still produced when the roots ji and f% coalesce. 

Finally we turn to /i, the largest of the three roots in figure 3. Here again, the root 
moves to the left as r decreases so that /i decreases as we move to the interior of the 
geometry indicating a positive mass. However, as noted above, this root coalesces with 
at r = ro and becomes complex for smaller values of r. Defining /i(ro) = /o and Taylor 
expanding eq. (4.19) about this point, we find 

/i(r)^/o + -^ ( r ^I±\ 1 ' 2 where 7 2 = -i/ i "(/o) = ^|/o-A. (4.20) 
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Figure 3: Graph of /i(/oo) for H = —0.005 and A = 0.145. The curve shifts down as r decreases. 
Here, the three curves correspond to w 4 /r 4 = 0, 1 and 1.752 from the top to bottom. The slope 
is negative for the roots, fi and fy, indicating that these are stable solutions while it is negative 
for /3 indicating the graviton is a ghost in this background. At w 4 /r 4 = 1.752, the roots /i and /3 
coalesce and a curvature singularity appears in both solutions. 



Further calculating the curvature using this result yields 

RabcdR abCd OC ^- - 1 (4.21) 

r L 4 (r - r^Y 

showing that the spacetime has a naked singularity at this point. One could again examine 
these solutions with cj 4 < 0. However, in this case, f± moves to the right, indicating a 
negative mass, and a naked singularity arises with fi — > +oo as r — > 0. 

This discussion shows that in case (g) from table 1, only the solution fi corresponds 
to an asymptotically AdS black hole. The other roots, f± and fs, are both asymptotically 
AdS but produce spacetimes with naked singularities. Examining the other cases in the 
table in a similar way, one finds that in each parameter regime with an AdS vacuum, there 
is a single black hole solution corresponding to the smallest positive root of eq. (4.19). The 
only exception is case (c) where we must also be in the regions denoted (co) or (ci). These 
restrictions are related to the possibility that a naked singularity will arise if the function 
h(f) is not monotonic in the range / G [0, foo], as explained for the root f\ in the example 
above - see also figure 4. First of all, in the region (co), p is negative and h(f) has no 
extrema at all. In regions (ci) and (C2), p > and so one must examine the extrema /o of 
h(f). In region (ci) to the left of (g) where < 0, both of the extrema /o > foo and so 
h(f) is monotonic in the desired range. On the other hand, one finds < / < in region 
(02) to the right of (g). Therefore in this parameter regime, the solution develops a naked 
singularity when r reaches the value where f = fo- However, the solution corresponds to 
a black hole with smooth event horizon for parameters in the regions (cq) and (ci). All of 
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Figure 4: The function h(f) plotted for /j, = —0.045 and A = 0.4, typical parameters in the region 
(02). The single (real) root f-y corresponds to a ghost-free asymptotically AdS solution. However, 
when the radius reaches a value where h'(f) = at the root, indicated with the black dot, the 
geometry becomes singular. 



our results with regards to which root yields a black hole solution are summarized in table 
1. 



4.2 Curved horizons 

As in eq. (2.16), we can again generalize the metric ansatz to include spherical and hyper- 
bolic, as well as planar, horizons: 

ds 2 = -(k+ ^f(r)) N(r) 2 dt 2 + ; - + r 2 dit (4.22) 

where d£ 2 is given by 



L 2 ') k + £f(r 



k = +l: dn% , 

k = : -3 (dx 2 + dy 2 + dz 2 ) , (4.23) 



k = -1 : cE? 



L 2 

J 3 ■ 



As in eq. (2.17), we have above the metric on a unit three-sphere for k = +1 and on a three- 
dimensional hyperbolic plane with unit curvature for k = — 1. The analysis at the beginning 
of section 4.1 follows through unchanged. Implicitly, we will assume N(r) 2 = 1/ foo but 
more importantly, the solutions are again determined by eq. (4.4): 



l-/ + A/ 2 + ^/ 3 = ^. (4.24) 
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Hence one arrives at the same solutions for / as previously found. The difference between 
planar and curved horizons is that the usual horizon equation g tt = now becomes / = 
—k^r- We have not made a complete analysis of the structure of the new spacetimes 
throughout fi-X plane but let us make the following preliminary remarks. 

We can develop a qualitative picture of the solutions using the same graphical approach 
as in the previous section. However, as well as following the behaviour of the h(f), we must 
now also keep track of the critical value of / where a horizon can form when k 7^ 0, i.e., 
fh = —kL 2 /r 2 . While the mass parameter oj 4 controls how quickly h(f) (and its roots) 
are shifting as r varies, the rate of change in fh is controlled by the cosmological constant 
scale L while the direction is controlled by k. Our first observation then is that if we have 
a black hole solution with k = 0, then for large masses oj S> L, we will always find the 
new (curved) horizon equation with ~ oj. That is, in this regime, the relevant root of 
h(f) moves much more quickly as r decreases than fh- Hence the root reaches / = at 
r = oj (as discussed for the planar horizons) while \fh\ = L 2 /uj 2 <C 1 and so there should 
be a nearby solution for the horizon condition. Therefore we expect that there are smooth 
black hole solutions with spherical or hyperbolic horizons in all of the same regions of the 

plane where they were found for planar horizons, in the previous section. However, in 
these regions, one may find that there is a lower bound on the mass of these black hole 
solutions different from oj 4 = 0. 

Let us consider small masses first for spherical horizons with k = +1 and f\ = —L 2 /r 2 . 
As r decreases, the smallest positive root of eq. (4.19) moves to the left, approaching / = 
as t — y oj. At the same time, fh starts at zero at r = 00 and then moves to the left 
to negative values as r decreases. Hence to form a horizon, the root must 'catch up' to 
fh- If we tune oj 4 to smaller and smaller values, slowing down the rate at which the root 
moves, it becomes clear that it may never coincide with fh- For example, consider cases 
(d) with p > 0, (e), (f) and (h) in table 1. In each of these cases, there will be a value 
fo < for which h'(fo) = 0, where the spacetime develops a singularity as described 
in section 4.1 - note that this previous discussion does not change for k 7^ 0. Hence if 
fh(oj) = —L 2 /oj 2 < /o, then the root will not be able to reach f\ before hitting /o. Hence 
in this situation, the spacetime will contain a naked singularity. In cases (c), (d) with 
p < and (g), h(f) is monotonic for negative values of /. However, now for small r (i.e., 
r <C oj, L) the root behaves as / ~ — [a; 4 / r 4 )] 1 / 3 . Hence the root is growing much more 
slowly than fh oc r~ 2 in this regime and again it becomes apparent that the root will never 
catch up to fh- Hence the solution will again have a naked singularity as r — > 0. Thus 
our final conclusion is that for spherical horizons with k = +1 there will always be a lower 
(positive) bound on the mass parameter below which no black hole solutions exist. This 
is, of course, qualitatively, the same result as found for spherical black holes for Einstein 
gravity with a negative cosmological constant [25] . We have not calculated the exact value 
of the lower bound but one must find oj ~ L, with the precise proportionality constant 
determined by the gravitational couplings, A and fi. 

We might add that since fh < with k = +1, the discussion given in the previous 
section remains unchanged for the solutions asymptotic to the ghosty vacua, the extra 
stable AdS vacuum in case (g) and the AdS vacuum in case (02). That is, there will be no 
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spherical black hole solutions in any of these cases. 

Now let us turn to considering small masses for hyperbolic horizons with k = — 1 and 
f h = L 2 /r 2 . In this case, fh moves to the right to positive values as as r decreases and so 
more exotic possibilities arise. For example, even if u = and the root does not move, fh 
will increase and eventually coincide with the root. That is, a horizon will appear even if 
the mass is set to zero and by continuity, we must also find black holes with small negative 
masses. Note that with w 4 , the root moves to the right to values larger than the initial foo. 
Here, cases (e), (f) and (g), as well as the region denoted (ci), are distinguished because 
h(f) is not monotonic for positive values of / > f^. Hence the negative mass solutions will 
only contain a horizon if fh catches up to the root before reaching the point where h! = 0. 
This will set a (negative) lower bound on the mass with w 4 ~ — L 4 where again the precise 
bound will be determined by A and fi. Cases (co), (d) and (h) with k = — 1 are even more 
striking. For these cases, h(f) is monotonic for / > f^ and so a singularity only develops 
as r — > and / — > oo. However, in this regime, the root grows as / ~ [|c<J 4 |/(|/i| r 4 )] 1 / 3 
while fh = L 2 /r 2 . Hence the solution will also contain a horizon with r 2 ~ I^IL 6 /^ 4 ) 
for arbitrary negative values of uj 4 . Hence we have found that in any of the parameter 
regimes where planar black holes exist, there will be black holes with hyperbolic horizons 
with negative masses. In cases (ci), (e), (f) and (g), there is a lower bound on how negative 
the mass can become but in cases (co), (d) and (h), the hyperbolic black holes can have 
an arbitrarily large negative mass. These results are not entirely surprising given that 
hyperbolic black holes with negative masses exist for for Einstein gravity with a negative 
cosmological constant [26], although there is a lower bound on the mass there. 

We might briefly also consider the effect of setting k = — 1 in the cases where no planar 
black hole solutions could be found, i.e., , the solutions asymptotic to the ghosty vacua, the 
extra stable AdS vacuum in case (g) and the AdS vacuum in case (02). It is clear that even 
in these cases a smooth horizon forms with u = since fh = L 2 /r 2 moves to positive values 
as r decreases. Again by continuity, black hole solutions will also exist for small positive 
and negative values of a; 4 . As above, the largest root will grow as / ~ [w 4 /(^r 4 )] 1 / 3 in 
a regime where r — > 0. Since the critical value fh grows more quickly, we conclude that 
the solutions which asymptote to the ghosty vacua in cases (e) and (f) (with jjl > 0) will 
form a horizon for arbitrarily large positive values of uj 4 . Similarly, the solutions which 
asymptote to the second stable AdS vacuum in case (g) or to the AdS vacuum in case 
(C2) will form hyperbolic horizons for arbitrarily large negative masses. Hence we find that 
exotic hyperbolic black holes also exist in parameter regimes even where no planar black 
holes formed. 

4.3 Higher dimensions 

With the construction described in section 3.1, we extended quasi-topological gravity to 
higher dimensions D > 7 - recall that our new curvature-cubed interaction did not affect 
the equations of motion in D = 6. The general action for D > 7 (as well as D = 5) is: 
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(4.25) 



where A4 and Z D are given in eqs. (2.1) and (3.10), respectively. We choose the metric 
ansatz in eq. (2.16), which can describe D-dimensional black holes with planar, spherical 
or hyperbolic horizons (for k = 0, +1 and -1, respectively). The coefficients of the action 
(4.25) are chosen so that substituting in this metric (2.16) yields 



/ = 



16vrG T 



J d° X {D l ] D N{r) [r^ (l - / + A/ 2 + »f)]' . (4.26) 



The lapse must again be constant and a convenient choice is N(r) 2 = I//00 as described in 
section 2 - with the implicit assumption that /oo > 0. If we wish to consider the vacuum 
solutions, f(r) is also fixed to be constant with f(r) = where these solutions are again 
determined by eq. (4.12). Hence the various vacuum solutions are again distributed as 
described in table 1. Assuming > 0, the D-dimensional metric becomes 

ds 2 = — (k + f- + — ^— + r 2 dt 2 (4.27) 

V L J loo k + -j^/oo 

where dt\ is given in eq. (2.17). These solutions correspond to a spherical (k = +1), flat 
(k = 0) or hyperbolic (k = —1) foliation of AdS D . From the g rr component, we see that 
the radius of curvature of the AdS D spacetime is 

L 2 = L 2 //oo • (4.28) 

In considering the black hole solutions, the only difference from the analysis for D = 5 
in the previous sections is that we replace: w 4 /r 4 — > oj ^ 1 jr D ~ x . In particular, the latter 
substitution is made in q in eq. (4.6). Further, eq. (4.19) is replaced with 

Hf) = (l - ^) ~f + \f 2 + vf = 0. (4.29) 

The remainder of the analysis and the results in section 4.1 carries over unaltered. In 
particular, table 1 correctly describes the planar black hole solutions for D > 7. Similarly, 
the discussion of black holes with curved horizons in section 4.2 remains largely unchanged. 
In certain cases, the discussion for small masses referred to the behaviour of the root, 
now, of eq. (4.29) as r approaches zero. In the present case, this behaviour changes to 
/ ~ [uj ^ 1 /(/ir ^ 1 )] 1 / 3 while the behaviour of the critical value remains fh = L 2 /r 2 . 
Hence the root is grows more slowly than fh for D = 5, at the same rate for D = 7 and 
more quickly for D > 8. As a result, one finds, in cases (c), (d) with p < and (g), that 
there is a positive lower bound for the mass of spherical black holes for D = 5 and 7 but the 
lower bound is simply uj ^ 1 = for D > 8. Further, the conclusion that hyperbolic black 
holes exits in cases (co), (d) and (h) with arbitrarily large negative masses only applies 
for D = 5. Similarly, some of the details about the formation of hyperbolic horizons, in 
regions of the coupling space where planar black holes do not exist, change depending the 
spacetime dimension D. 
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As noted above, this discussion applies for D > 7, but in those dimensions, one already 
has a cubic order theory in Lovelock gravity, which also reproduces eq. (4.29). Hence for 
D > 7, the black holes for quasi-topological gravity discussed here would be the same as 
those in cubic Lovelock theory, as considered in [8, 27]. 



5. Black Hole Thermodynamics 

We now turn to the thermodynamic properties of the black hole solutions of quasi-topological 
gravity. Our focus will be on the planar (k = 0) black holes in five dimensions. The exten- 
sion of these results to curved horizons and higher dimensions is straightforward. 

First, we use the standard approach to calculate temperature: analytically continue 
the metric to Euclidean signature with t = —it and periodically identify r to produce an 
everywhere smooth Euclidean section. Interpreting the period of r as the inverse temper- 
ature, we find 

rp 1 r hf'\rh 



4vr L 2 ^f 



(5.1) 



which assumes the lapse is chosen as N = 1/V/oo- For the planar black holes, = u and 
further we will use f(r = co) = 0. Then we can evaluate the f'\ rh by differentiating the 
constraint equation (4.4) and evaluating the result at r = co. A simple calculation yields 
/'!« = £, giving: 



(5.2) 



In passing we note that extending this calculation to general dimensions and curved 
horizons yields 
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(5.3) 



where the precise location of the horizon must still be determined for a spherical or 
hyperbolic horizon but, of course, = uj in the planar case (k = 0). 

Next, we calculate the entropy and energy densities of the black holes following the 
Euclidean action approach, as already sketched for GB theory in section 2. That is, we 
identify the Euclidean action for the black hole solution, as the leading contribution to the 
free energy, i.e., Ie — F/T. Evaluating the Euclidean action yields 
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Here V3 = J d 3 x is the regulator volume for the (spatial) gauge theory directions and as 
usual, we have set TV" = l/V/oo- The final result is simplified using the constraint (4.4) to 
produce the following asymptotic expansion of f(r): 

uj 4 1 

/ " /o0 "^(l-2A/ 00 -3 / x/^) + ---- (5 - 5) 

Keeping only the divergent and finite terms in the limit Rq — > 00, eq. (5.4) then reduces 
to: 



Ie[T] 
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We remove the divergence in this expression by subtracting the action for the AdS vacuum 



where T' is chosen so that the periodicity of the AdS background matches that of the black 
hole at the regulator surface r = Rq: 



T' T vToo T V 2R 4 f 00 (l-2\f 00 -3fifl)J ■ 
Combining these expressions yields a simple expression for the free energy: 

m = T [I E [T] - I% [T] ) = -nggjjr; = 



(5.8) 



(5.9) 



where we have implicitly taken the limit Rq — > 00 above. As noted in section 2, for planar 
AdS black holes, we do not have to account for the generalized Gibbons-Hawking term 
[20] or the boundary counter-terms [21] in this calculation of the free energy. Finally, to 
eliminate the regulator volume, we work with the free energy density: T = F/V3. Then 
we calculate the energy and entropy densities as 

^ 2 d 3(tt^L) 3 T 4 



P--T^/T)= - — / , (5.10) 
dF _ (vry^LT) 3 _ 1 W 3 

Note that these expressions are 'identical' to those appearing in eq. (2.14) for GB gravity, 
however, implicitly depends on the additional coupling fi. We also comment that 
these results obey the relation p = |Ts, as expected for a four-dimensional CFT at finite 
temperature. 

5.1 Noether Charge Approach to Entropy Density 

In this section, we verify that the entropy density in eq. (5.11) matches the Wald entropy 
[22]. Of course, we find agreement [28] but the results for the Wald entropy are also readily 
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extended to higher dimensions and curved horizons. Wald's prescription for the black hole 
entropy is 

S = -2tt (f d D - 2 xVhY abcd e ab e cd where Y abcd = -|^- (5.12) 

J uRabcd 

C is the Lagrangian and i a t, is the binormal to the horizon. For the static black holes 
considered here, Y = Y abcd i a j J i a j j is constant on the horizon and so the entropy is given 
simply as 

S = -2ttYA, (5.13) 

where A = <f> dP~ 2 x\fh is the 'area' of the horizon. Let us divide up the terms in the 
action (4.25) according to their powers of the curvature tensor: the Einstein term, the GB 
interaction and the quasi-topological interaction. (Of course, one also has the cosmological 
constant term but it does not contribute in Wald's formula (5.12).) Following the above 
prescription, as usual, the Einstein term yields Y\ = — 1/(8ttG d ) and the resulting entropy 
is the expected Bekenstein-Hawking entropy S = A/(4G D ). Applying this formalism to 
GB terms, we find 

* ~ -lk (D-WO-4) < R " 2 + ^ + 2 *""> • (5 ' I4) 

where this expression applies for a general static black hole metric. One can think that 
the tensor components above are presented in an orthonormal frame or alternatively in a 
coordinate frame, as long as the indices are in precisely the raised and lowered positions as 
shown. Integrating this over the horizon gives the following contribution to the entropy as 

A XL 2 
S2 = 2G£(£>-3)(£>-4) (jR " 2 ^ + K) + 

'' 2 I 2A./V;,i;. . (5.15) 



4G D D - 4 



where we are using the general metric (2.16) with iV 2 = l//<x> to evaluate the second line. 
Finally turning to the quasi-topological contribution in the action, we find 



1 8(2D-3)/iL 4 

Y3 - 



1 ( Titm r>rn _ r>tm nr n\ /r ic\ 

4itG d (D -6)(D — 3) (3D 2 - 15D + 16) |_ 2 [K tnK rm " ™K mt ) l^.lbj 

_i_q„ Tftrmn p i „ [ pir p rn _ ptr p rn , \ ( p jymnpr , p T>mnpt\ \ 

-t-ot-2-n. J^trmn T1.3 I Ji tm-^r 11 rm^t 1 \ jr >mnpr rL t f^mnpt 1 ^ j I 

[ OR E?' r _l ^ J? Tjmnpq] . c 5 I r>t pr pt pr , pr r>mn , pt r>mn\ 

+C4 I Zil it tr + -Umnpglt I + — ^it til r — it r it t + it mrn K + it m t ra il ) 

+ ^C 6 (R rm R rm + R tm Rtm) + (RmnR mn + R {R T r + R*t)) + T^ 8 ^ 

We have left the coefficients arbitrary above but it is understood that they are to be fixed 
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as in eq. (3.9). Now integrating over the horizon yields 



A 8(2L>-3)//L 4 

^3 =~ 



1 / Tjtm r>rn _ pirn pr n\ 
2 x 1 ^ tn 1 ^ rm lx rn 1 *- mt ) 



2G D (D -6)(D- 3) (3D 2 - 15D + 16) 

J-1r> Tjtrrnn p , [ ptr p m _ ptr T} m \}_ ( J} nrnnpr , p jjmnpt\ \ 

-TOC2-11 JM.rmn > c 3 I - rl tm^r 1 ' rm r ^t T ^ \ i - f lmnpr- a ~r J^mnpt 1 ^ j I 

I OR f?* r I ^ p nmnpq \ . c 5 f pi pr pt pr . pr pmn , pt pmn\ 

3 C 7 3 

+ ^C 6 (R Tm R rm + i? m Rtrn) + — (RmnR mn + -R (^r + -R t)) + T^ 8 "^ 

^^f^/W). < 5 ' 17 > 

Combining all of these expressions, we arrive at the Wald entropy for quasi-topological 
gravity: 

A ( 2(D-2) 3(D-2) 2 



Evaluating this expression on a planar horizon yields the simple result, S = A/ (4G D ), i.e., , 
the higher curvature contributions vanish on planar horizons. If we divide by the regulator 
volume, this yields the entropy density: 

which, of course, agrees with the result in eq. (5.11) for D = 5. For the case of black holes 
with curved horizons, as in section 4.3, we find that the entropy is given by 

A / 2(D-2) L 2 3(D-2) l2 L 4 \ . . 



6. Equations of Motion 

We have found that the equations of motion for quasi-topological gravity take an incredibly 
simple form with the ansatz (2.16) for a static AdS black hole - e.g., see eq. (4.26). In this 
section, we would like to investigate the equations of motion in greater generality to gain 
some further insight into this simplicity. In particular, we will see the linearized equations 
of motion of graviton fluctuations also exhibit a certain simplicity. We have already argued 
that the new cubic interactions constructed in section 3 do not have a topological origin - 
see appendix A. Hence this cannot be the source of the simplicity noted above. 

Let us begin with the cubic-curvature interactions in eq. (3.1). First we set c u = = c 13 
and then find the general contribution these terms would make to the metric equations of 
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motion: 
_L 5J_ 

,ab 



c 8 (-3R 2 R ab + 3 (R 2 ) ;ab ) + c 7 ( — R c R d C R ab — 2RRacR C b + (Rc Rd, C );ab ~ D (RRab) 
+2{RR c [a ). b)c ) + c 6 (-3R ac R c d R d b - ^U(R ac R c b ) + 3(R c d R d (b ) . a)< \ 

+C 5 ^-3i? c ( (I i?5) rf c e i?i?' ie + 2(i?( a | c d e ii ce )| - □(i? ac 5 (i i? cd ) + (i?( a C R b ) d )- c d 

-(R cd R a b):cdj + c 4 (^-2RR acde R b cde - R ab R cde fR cde f + {R cde fR cde ^);ab 

+^{RRa Cd b);cdj + C 3 (^-2R acd e R b cd f R e f - R de f C R de f( a R b y + (R cde f R Cde (a);b)f 
+ -^{R Cde aRcdeb) + 2(R edC ( a Rb)e);cd + ^{R C {ab) C Re d ) -cd^j + C 2 ^Rac^ Rde^ 9 R f 

— 3R c def R C9£ aRj gb d + 3(R e d f( b R a ) eC ^);cd 

-3(R (a e b /R e d f c ), dc ) + g ab (c, Qi? 3 - 3D(i? 2 )) + c 7 Qi?i? c d i?/ - (i?i? cd ) ;c(i 
-□(#/#/))+ c 6 Qi?/i?/i? e c - ^(i? c eJ R ed ) ;c ^ +c 5 (±R cd R cedf R e f 
-(R e c f d R e f).Jj+c 4 (^RR cdef R cde f - a(R cdef R cde ^ + c 3 {^R c ^R cd JR fg 

-^(fl^Ww) + \c 2 R cd ^R e f h R gh cd + ^AViV/i^v) . (6.1) 

Since we have eliminated the terms involving derivatives of the curvature from the action by 
setting c u = = c 13 , the above contributions contain at most terms with four derivatives 
of the metric, such as in □ {RR ab ). This may be misleading, however, since we have made 
no attempt to simplify this expression using, e.g., the Bianchi identities. A useful check is 
to choose the coefficients above so that eq. (3.1) corresponds to the six-dimensional Euler 
density (3.7), i.e., choose c x = — 2c 2 , c 2 = 4 and the remaining coefficients as in eq. (3.3). 
In this case, we were able to verify that any terms involving derivatives of curvatures can 
be eliminated and the expected field equations involving only factors of the curvature were 
reproduced - e.g., see [27, 29]. Indeed Lovelock's general discussion [14] dictates that Xq is 
the only gravity Lagrangian which is cubic in curvatures for which the equations of motion 
do not contain any derivatives of curvatures. Hence one must expect equations of motion 
that include terms with derivatives of curvatures for any other choice of the coefficients 
q and, in particular, for quasi-topological gravity with C; as in eq. (3.9). While we still 
made some effort to 'tidy up' the derivative terms in eq. (6.1), we did not produce any 
particularly illuminating results. 

As a next step, we examine the linearized equations of motion for a graviton pertur- 
bation in quasi-topological gravity. Hence we fixed the coefficients as in eq. (3.9) and then 
substitute into the above expression (6.1): g ab = g^ + h ab where g [ ° ] b is a solution of the 
full equations of motion. Again, the resulting expression of a generic fluctuation is rather 
complicated and so to proceed further, we restrict our attention to transverse traceless 
gauge with V a h ab = 0) and h a a = 0. This choice simplifies the result somewhat and the 
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four-derivative contribution is proportional to: 

(D - A)R cdef h de . cf(ab) + R cd (Dh cd ). (ab) - 2R cd {Uh c{a ). b)d + ^^R cd iUh ab ), cd 
+2 (Uh\ a ) ' de R lcde \ b) + g ab (Dh cd ). ef R ced f + 2R c[a U 2 h c b) - X -g ab R cd U 2 h cd 
+ (D - 3) n 2 h cd R cadb - 1 -J^- Y) U 2 h ab , (6.2) 

where we use the standard notation T( a b) = \ {T ab + T ba ). Hence despite the reduction 
in the number of terms, we see that for general backgrounds, the linearized equations of 
motion for (physical) gravitons include four-derivative contributions. 

Of course, if we were considering gravitons propagating in flat space, these contribu- 
tions would all vanish because R a bcd = in the background spacetime. The same result 
would apply for any interactions which are cubic in curvatures for a flat background. How- 
ever, in the present context, it is natural to consider gravitons propagating in an AdS^i 
background. In this case, of course, the background curvature is nonvanishing and so 
one might expect these terms (6.2) will still appear in the linearized equations of motion. 
However, further simplifications can be expected since AdS^) is a maximally symmetric 
spacetime with 1 

Rabcd = — -~r (g ac gbd - g a dgbc) ■ (6.3) 

Remarkably, one finds that upon substituting this background curvature into eq. (6.2), all 
of the remaining four-derivative terms cancel! Further, because the background curvature 
(6.3) is covariantly constant, there are no nontrivial terms with only three derivatives acting 
on the graviton. Therefore, in quasi-topological gravity, the linearized graviton equation 
in an AdS/j background is only a second-order equation. 

With this result in hand, we next construct the linearized equation of motion for the 
graviton in the AdS/j vacuum solutions. Hence we must consider the full action (4.25) for 
quasi-topological gravity, including the cosmological constant, Einstein and Gauss-Bonnet 
terms as well. The equation of motion for the transverse traceless graviton can then be 
written as: 

87rG D f ab . (6.4) 



-l(l-2A/ 00 -3/ i /^) 



V2 1 i ^Joc i 
Kb + -pr h ab 



L 2 

We have added a stress-tensor on the right-hand side, as might arise from minimally cou- 
pling the metric to additional matter fields or from quadratic or higher order contributions 
in the graviton. The second bracketed factor on the right-hand side is the standard Ein- 
stein equation of motion for gravitons in an AdS background with curvature L = L/y/ 
[30]. One can recognize the first bracketed factor as the slope of the cubic equation (4.12) 
determining f^, i.e., see eq. (4.18). Hence this slope determines the sign of the graviton 
propagator or alternatively the sign of the coupling of the graviton to the stress tensor. 
We see that the appropriate sign for a well-behaved graviton is negative since this factor 



1 Note that here we are distinguishing L, the curvature scale of the AdS background, from L, the AdS 
length scale appearing in the action. In particular, recall that in the vacuum solutions above, we found 
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reduces to -1 when A = = \i. Hence as discussed in section 4.1, the AdS vacua are only 
stable when the slope is negative while the graviton is a ghost in AdS backgrounds where 
the slope is positive. 

While restricting to transverse traceless gauge is a convenient simplification to deter- 
mine the linearized equations in an AdS background, we can do better. A straightforward 
argument shows that the above gauge-fixed equations of motion (6.4) extend to the full 
linearized Einstein equations [30] : 2 



where g } is the background AdS metric. We know that the full linearized equations come 
from a covariant expression and hence they must invariant under the 'gauge' transforma- 
tions: Sh a b = V a £b + Vfe£ a . Now the linearized Einstein equations (6.5) are certainly 
invariant under these transformations and reduce to eq. (6.4) upon fixing to transverse 
traceless gauge. However, there may be additional contributions which are both gauge 
invariant and completely vanish for transverse traceless modes. For example, the full equa- 
tions may include an additional contribution proportional to 



which is gauge invariant but would not contribute in eq. (6.4). In this particular case, it 
is easy to argue that such a contribution could not arise from the variation of an action 
(quadratic in h a b). However, given the action (4.25) for quasi-topological gravity, another 
approach is to evaluate the quadratic action (using Mathematica) and subsequently exam- 
ining the equations of motion for some specific trial perturbations which are not transverse 
or traceless. In every case considered, the latter equations match precisely the results 
expected from eq. (6.5). This clearly shows that there are no additional terms of the 
form given in eq. (6.6) but more importantly that there are no additional four-derivative 
contributions in the full linearized equations without gauge-fixing. The fact that gravi- 
tons propagating in an AdS background simply obey the same equations of motion as in 
Einstein gravity plays an important role in understanding the holographic properties of 
quasi-topological gravity [10]. 

7. Discussion 

In section 3, we have constructed a new gravitational action which includes terms cubic 
in the curvature. Our construction was motivated by the simple equations (2.5) arising 
to determine the black hole solutions in GB gravity. We were able to reproduce a similar 
structure (4.4) for our new theory. We wish to emphasize how remarkable this result is. 

2 We would like to thank Miguel Paulos for discussions on this point and confirming that several trial 
perturbations satisfied eq. (6.5) with Mathematica. 
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In section 6, we showed that the full equations are fourth order in derivatives. However, 
in section 4, we found that once the geometry of the horizon is fixed, the static black hole 
solutions are fixed by a single integration constant! It seems that the symmetry imposed on 
the background geometry must play an important role in producing the simplicity of these 
solutions. While our approach was to substitute the ansatz into the action, it would be 
interesting to work directly with the equations of motion and formalize this result in terms 
of a 'Birkhoff theorem.' It would also be interesting to see to what extent this simplicity 
extends to spinning or electrically charged black holes in quasi-topological gravity. 

We also saw that the linearized equations of motion for gravitons propagating in the 
AdS backgrounds reduced to the same second order equations as for Einstein gravity in 
section 6. There we saw more or less directly that this simplification comes about due to 
the maximal symmetry of the AdS spacetime. On the other hand, we should not expect 
that the four-derivative contributions (6.2) to the equations of motion cancel in the black 
hole backgrounds. Hence the quasinormal spectrum of the black hole solutions should be 
studied in detail. In particular, this spectrum may reveal that these solutions are unstable 
for certain values of the gravitational couplings. 

The simplicity of the graviton equations in an AdS background has the interesting 
consequence that the standard holographic rules apply in matching the metric fluctuations 
to the stress tensor of the CFT. In a general higher derivative theory, implicitly the gravi- 
ton would be matched with some higher dimension operator, as well as the stress tensor 
[6]. This interpretation arises because the higher derivative equations allow the metric 
fluctuations to have more that the standard asymptotic behaviour in the AdS geometry. 
In any event, this complication is evaded in quasi-topological gravity and so the effect of 
the higher derivative gravitational terms will only be felt in the higher n-point couplings 
of the stress tensor. 

Having motivated the construction of quasi-topological gravity by considerations of 
the AdS/CFT correspondence, one might ask how the universality class of dual CFT's has 
been expanded. The gravitational theory is defined by three independent dimensionless 
parameters: A, fi and L D ~ 2 /G D . In general, the three-point function of the stress tensor of a 
CFT in four or higher dimensions is also characterized by three independent (dimensionless) 
parameters [31]. This match in the counting of these parameters is not a coincidence, as the 
discussion of [32] that holographically modelling the full range of these CFT parameters 
requires the introduction of curvature-squared and curvature-cubed interactions in the bulk 
gravity theory. We make precise the mapping between the gravitational couplings and the 
dual CFT parameters, as well as exploring other holographic aspects of quasi-topological 
gravity in [10, 33]. 

Another natural extension of this work is to consider analogous gravitational interac- 
tions with higher powers of the curvature. Of course, Lovelock gravity provides an infinite 
sequence of (curvature)™ interactions which still allow for a certain calculational control 
in the gravitational theory. However, because of their topological origin, these Lovelock 
terms only contribute to the equations of motion for D > 2n + 1, i.e., in the context of 
the AdS/CFT correspondence, for dual CFT's with d > In. However, our expectation is 
that our construction of five-dimensional quasi-topological gravity with curvature-cubed 
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interactions can be extended general (curvature)™ interactions and this has been verified 
by some preliminary calculations [9]. In fact, our conjecture is that in D > 2n + 1, an in- 
dependent interaction can be constructed for each independent scalar contraction of (Weyl 
tensor)™. For D < 2n + 1, Schouten identities will reduce the number of independent inter- 
actions, as seen in the present analysis of the curvature-cubed interactions. However, our 
understanding of this issue remains incomplete, as we are still uncertain as to why there 
was no effective curvature-cubed interaction in six dimensions. In any event, better under- 
standing the number of independent quasi-top ological gravity terms with higher powers of 
the curvature will be an interesting direction of study. Important new insights into this 
question were given in refs. [12, 13]. There it was shown that the curvature-cubed interac- 
tions constructed here can be simply expressed in terms of scalar contractions of the Weyl 
tensor combined with the six-dimensional Euler density. Their construction immediately 
generalizes to an infinite family of higher curvature interactions with similar properties. 
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A. A New Topological Invariant? 

In section 3.1, our construction produced an nontrivial interaction (3.10) for any number 
of dimensions D > 5. However, we noted Zq, the six-dimensional expression, did not con- 
tribute to the equations of motion for the black hole metric (2.3) (extended to D = 6). In 
particular, this also means that the the AdS vacua are unaffected by the addition of Z§ 
to the gravitational action. This behaviour is reminiscent of Lovelock gravity where, for 
example, the Euler density Xq provides a nontrivial gravitational interaction for D > 7 but 
because of its topological origin, it leaves the equations of motion unaffected for D < 6. 
Hence one might be tempted to believe that Zq yields another topological invariant in six 
dimensions. However, in the following, we demonstrate that J cfix^/g Z§ is not a topolog- 
ical invariant by explicitly evaluating this expression for certain specific six-dimensional 
geometries. 

As our first test, we evaluate this expression on a deformed six-sphere with metric: 

ds 2 = R 2 [d9 2 + sin 2 6 (1 + a sin 2 6) n dftl] (A.l) 
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where n (implicitly an integer) and a are constants denning the deformation away from the 
round six-sphere. We then find: 



where we have normalized Xq as in eq. (3.7). Hence we see that both of these results are 
independent of the deformation parameters. Of course, for Xq, this occurs because the 
integrated expression is a topological invariant. While again suggestive for Z e , this result 
is by no means conclusive and hence we consider a further test. 

Next we consider the following metric in which the spheres in the direct product S 2 x S 4 
are deformed: 



where the deformation is characterized by the constants a and b. In this case, we find: 



where F(a,b, R/L) is a complicated (and not particularly illuminating) function of both 
deformation parameters and the relative radius of curvature of the two spheres. Hence this 
result makes clear that Zq does not yield a topological invariant. 

References 

[1] P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting quantum field 
theories from black hole physics," Phys. Rev. Lett. 94, 111601 (2005) [arXiv:hep-th/0405231]; 
P. Kovtun, D. T. Son and A. O. Starinets, "Holography and hydrodynamics: Diffusion on 
stretched horizons," JHEP 0310, 064 (2003) [arXiv:hep-th/0309213]. 

[2] A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," Phys. Rev. 
Lett. 93, 090602 (2004) [arXiv:hep-th/0311175]; 

A. Buchel, "On universality of stress-energy tensor correlation functions in supergravity," 
Phys. Lett. B 609, 392 (2005) [arXiv:hep-th/0408095]; 

P. Benincasa, A. Buchel and R. Naryshkin, "The shear viscosity of gauge theory plasma with 
chemical potentials," Phys. Lett. B 645, 309 (2007) [arXiv:hep-th/0610145]; 
D. Mateos, R. C. Myers and R. M. Thomson, "Holographic viscosity of fundamental matter," 
Phys. Rev. Lett. 98, 101601 (2007) [arXiv:hep-th/0610184]; 

K. Landsteiner and J. Mas, "The shear viscosity of the non-commutative plasma," JHEP 
0707, 088 (2007) [arXiv:0706.0411 [hep-th]]; 

N. Iqbal and H. Liu, "Universality of the hydrodynamic limit in AdS/CFT and the 
membrane paradigm," Phys. Rev. D 79, 025023 (2009) [arXiv:0809.3808 [hep-th]]. 

[3] For example, see: 

A. Buchel, J. T. Liu and A. O. Starinets, "Coupling constant dependence of the shear 
viscosity in N=4 supersymmetric Yang-Mills theory," Nucl. Phys. B 707, 56 (2005) 
[arXiv:hep-th/0406264]; 




(A.2) 



ds 2 = R 2 \d6 2 + sin 2 6 (l + a sin 2 ef d<p 2 



(A.3) 





(A.4) 



- 30 - 



P. Benincasa and A. Buchel, "Transport properties of N = 4 supersymmetric Yang-Mills 

theory at finite coupling," JHEP 0601, 103 (2006) [arXiv:hep-th/0510041]; 

A. Buchel, "Shear viscosity of boost invariant plasma at finite coupling," Nucl. Phys. B 802, 

281 (2008) [arXiv:0801.4421 [hep-th]]; "Resolving disagreement for r\js in a CFT plasma at 

finite coupling," Nucl. Phys. B 803, 166 (2008) [arXiv:0805.2683 [hep-th]]; 

R. C. Myers, M. F. Paulos and A. Sinha, "Quantum corrections to r]/s," Phys. Rev. D 79, 

041901 (2009) [arXiv:0806.2156 [hep-th]]. 

[4] Y. Kats and P. Petrov, "Effect of curvature squared corrections in AdS on the viscosity of the 
dual gauge theory," JHEP 0901, 044 (2009) [arXiv:0712.0743 [hep-th]]; 
A. Buchel, R. C. Myers and A. Sinha, "Beyond rj/s = 1/4tt," JHEP 0903, 084 (2009) 
[arXiv:0812.2521 [hep-th]]; 

R. C. Myers, M. F. Paulos and A. Sinha, "Holographic Hydrodynamics with a Chemical 
Potential," JHEP 0906, 006 (2009) [arXiv:0903.2834 [hep-th]]. 

[5] M. Brigante, H. Liu, R. C. Myers, S. Shenker and S. Yaida, "Viscosity Bound Violation in 
Higher Derivative Gravity," Phys. Rev. D 77 (2008) 126006 [arXiv:htp-th/0712.0805]; "The 
Viscosity Bound and Causality Violation," Phys. Rev. Lett. 100, 191601 (2008) 
[arXiv:0802.3318 [hep-th]]; 

A. Buchel and R. C. Myers, "Causality of Holographic Hydrodynamics," JHEP 0908, 016 
(2009) [arXiv:0906.2922 [hep-th]]. 

[6] D. M. Hofman, "Higher Derivative Gravity, Causality and Positivity of Energy in a UV 
complete QFT," Nucl. Phys. B 823, 174 (2009) [arXiv:0907.1625 [hep-th]]. 

[7] X. H. Ge and S. J. Sin, "Shear viscosity, instability and the upper bound of the 

Gauss-Bonnet coupling constant," JHEP 0905, 051 (2009) [arXiv:0903.2527 [hep-th]]; 

R. G. Cai, Z. Y. Nie and Y. W. Sun, "Shear Viscosity from Effective Couplings of 

Gravitons," Phys. Rev. D 78, 126007 (2008) [arXiv:0811.1665 [hep-th]]; 

R. G. Cai, Z. Y. Nie, N. Ohta and Y. W. Sun, "Shear Viscosity from Gauss-Bonnet Gravity 

with a Dilaton Coupling," Phys. Rev. D 79, 066004 (2009) [arXiv:0901.1421 [hep-th]]; 

J. dc Boer, M. Kulaxizi and A. Parnachev, "AdSt/CFT^, Gauss-Bonnet Gravity, and 

Viscosity Bound," arXiv:0910.5347 [hep-th]; 

X. O. Camanho and J. D. Edelstein, "Causality constraints in AdS/CFT from conformal 
collider physics and Gauss-Bonnet gravity," arXiv:0911.3160 [hep-th]; 

A. Buchel, J. Escobedo, R. C. Myers, M. F. Paulos, A. Sinha and M. Smolkin, "Holographic 
GB gravity in arbitrary dimensions," arXiv:091 1.4257 [hep-th]. 

[8] J. de Boer, M. Kulaxizi and A. Parnachev, "Holographic Lovelock Gravities and Black 
Holes," arXiv:0912.1877 [hep-th]; 

X. O. Camanho and J. D. Edelstein, "Causality in AdS/CFT and Lovelock theory," 
arXiv:0912.1944 [hep-th]. 

[9] Miguel Paulos, unpublished. 

[10] R. C. Myers, M. F. Paulos and A. Sinha, "Holographic studies of quasi-topological gravity," 
arXiv: 1004.2055 [hep-th]. 

[11] A. Sinha, "On the new massive gravity and AdS/CFT," arXiv:1003.0683 [hep-th]. 

[12] J. Oliva and S. Ray, "A new cubic theory of gravity in five dimensions: Black hole, Birkhoff 's 
theorem and C- function," arXiv: 1003.4773 [gr-qc]. 



- 31 - 



[13] J. Oliva and S. Ray, "A Classification of Six Derivative Lagrangians of Gravity and Static 
Spherically Symmetric Solutions," arXiv: 1004.0737 [gr-qc]. 

[14] D. Lovelock, "The Einstein tensor and its generalizations," J. Math. Phys. 12, 498 (1971); 
Acquationes Math. 4, 127 (1970). 

[15] B. Zwiebach, "Curvature Squared Terms And String Theories," Phys. Lett. B 156 (1985) 315. 

[16] D. G. Boulware and S. Deser, "String Generated Gravity Models," Phys. Rev. Lett. 55 
(1985) 2656. 

[17] J. T. Wheeler, "Symmetric Solutions To The Gauss-Bonnet Extended Einstein Equations," 
Nucl. Phys. B 268 (1986) 737; 

J. T. Wheeler, "Symmetric Solutions To The Maximally Gauss-Bonnet Extended Einstein 
Equations," Nucl. Phys. B 273 (1986) 732; 

R. C. Myers and J. Z. Simon, "Black Hole Thermodynamics in Lovelock Gravity," Phys. Rev. 
D 38 (1988) 2434; 

R. C. Myers and J. Z. Simon, "Black Hole Evaporation and Higher Derivative Gravity," Gen. 
Rel. Grav. 21, 761 (1989). 

[18] R. G. Cai, "Gauss-Bonnet black holes in AdS spaces," Phys. Rev. D 65 (2002) 084014 
[arXiv:hep-th/0109133]. 

[19] S. Nojiri and S. D. Odintsov, "Anti-dc Sitter black hole thermodynamics in higher derivative 
gravity and new conhning-deconfining phases in dual CFT," Phys. Lett. B 521 (2001) 87 
[Erratum-ibid. B 542 (2002) 301] [arXiv:hep-th/0109122]; 

Y. M. Cho and I. P. Neupane, "Anti-de Sitter black holes, thermal phase transition and 
holography in higher curvature gravity," Phys. Rev. D 66 (2002) 024044 
[arXiv:hep-th/0202140]; 

I. P. Neupane, "Black hole entropy in string-generated gravity models," Phys. Rev. D 67 
(2003) 061501 [arXiv:hep-th/0212092]; 

I. P. Neupane, "Thermodynamic and gravitational instability on hyperbolic spaces," Phys. 
Rev. D 69 (2004) 084011 [arXiv:hep-th/0302132]. 

[20] R. C. Myers, "Higher Derivative Gravity, Surface Terms and String Theory," Phys. Rev. D 
36, 392 (1987). 

[21] V. Balasubramanian and P. Kraus, "A stress tensor for anti-de Sitter gravity," Commun. 
Math. Phys. 208, 413 (1999) [arXiv:hep-th/9902121]; 

R. Emparan, C. V. Johnson and R. C. Myers, "Surface terms as counterterms in the 
AdS/CFT correspondence," Phys. Rev. D 60, 104001 (1999) [arXiv:hcp-th/9903238]; 
R. B. Mann, "Misner string entropy," Phys. Rev. D 60, 104047 (1999) 
[arXiv:hep-th/9903229]. 

[22] R. M. Wald, "Black hole entropy is the Noether charge," Phys. Rev. D 48, 3427 (1993) 
[arXiv:gr-qc/9307038]; 

V. Iyer and R. M. Wald, "Some properties of Noether charge and a proposal for dynamical 
black hole entropy," Phys. Rev. D 50, 846 (1994) [arXiv:gr-qc/9403028]; 
T. Jacobson, G. Kang and R. C. Myers, "On Black Hole Entropy," Phys. Rev. D 49, 6587 
(1994) [arXiv:gr-qc/9312023]. 

[23] T. Eguchi, P. B. Gilkey and A. J. Hanson, "Gravitation, Gauge Theories And Differential 
Geometry," Phys. Rept. 66, 213 (1980). 



- 32 - 



[24] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large N field theories, 
string theory and gravity," Phys. Rept. 323, 183 (2000) [arXiv:hep-th/9905111]. 

[25] S. W. Hawking and D. N. Page, "Thermodynamics Of Black Holes In Anti-De Sitter Space," 
Commun. Math. Phys. 87, 577 (1983); 

E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2, 505 (1998) [arXiv:hep-th/9803131]. 

[26] R. Emparan, "AdS / CFT duals of topological black holes and the entropy of zero-energy 
states," JHEP 9906, 036 (1999) [arXiv:hep-th/9906040]. 

[27] M. H. Dehghani and R. B. Mann, "Thermodynamics of rotating charged black branes in 
third order Lovelock gravity and the counterterm method," Phys. Rev. D 73, 104003 (2006) 
[arXiv:hep-th/0602243]; 

M. H. Dehghani and R. B. Mann, "Thermodynamics of rotating charged black branes in 
third order Lovelock gravity and the counterterm method," Phys. Rev. D 73, 104003 (2006) 
[arXiv:hep-th/0602243] ; 

M. H. Dehghani and R. Pourhasan, "Thermodynamic instability of black holes of third order 
Lovelock gravity," Phys. Rev. D 79, 064015 (2009) [arXiv:0903.4260 [gr-qc]]; 
M. H. Dehghani and M. Shamirzaie, "Thermodynamics of asymptotic fiat charged black holes 
in third order Lovelock gravity," Phys. Rev. D 72, 124015 (2005) [arXiv:hep-th/0506227]. 
S. H. Hcndi and M. H. Dehghani, "Taub-NUT Black Holes in Third order Lovelock Gravity," 
Phys. Lett. B 666, 116 (2008) [arXiv:0802.1813 [hep-th]]. 

[28] V. Iyer and R. M. Wald, "A Comparison of Noether charge and Euclidean methods for 
computing the entropy of stationary black holes," Phys. Rev. D 52, 4430 (1995) 
[arXiv:gr-qc/9503052]. 

[29] Y. Decanini and A. Folacci "Irreducible Forms for the Metric Variations of the Action Terms 
of Sixth-Order Gravity and Approximated Stress-Energy Tensor," Class. Quant. Grav. 24 
(2007) 4777 [arXiv:hcp-th/0706.0691] 

[30] For example, see: 

H. Liu and A. A. Tseytlin, "D = 4 super Yang-Mills, D = 5 gauged supergravity, and D = 4 
conformal supergravity," Nucl. Phys. B 533, 88 (1998) [arXiv:hep-th/9804083]; 
G. Arutyunov and S. Frolov, "Three-point Green function of the stress-energy tensor in the 
AdS/CFT correspondence," Phys. Rev. D 60, 026004 (1999) [arXiv:hep-th/9901121]. 

[31] H. Osborn and A. C. Petkou, "Implications of Conformal Invariance in Field Theories for 
General Dimensions," Annals Phys. 231, 311 (1994) [arXiv:hep-th/9307010]; 
J. Erdmenger and H. Osborn, "Conserved currents and the energy-momentum tensor in 
conformally invariant theories for general dimensions," Nucl. Phys. B 483, 431 (1997) 
[arXiv:hep-th/9605009]. 

[32] D. M. Hofman and J. Maldacena, "Conformal collider physics: Energy and charge 
correlations," JHEP 0805, 012 (2008) [arXiv:0803.1467 [hep-th]]. 

[33] R. C. Myers and A. Sinha, "Seeing a c-theorem with holography," arXiv:1006.1263 [hep-th]. 



- 33 - 



